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This paper presents the study of the flow structure and heat transfer, and also their correlations on the four walls 


of a radial cooling passage model of a gas turbine blade. The investigations focus on heat transfer and aerody- 


namic measurements in the channel, which is an accurate representation of the configuration used in aeroengines. 


Correlations for the heat transfer coefficient and the pressure drop used in the design of radial cooling passages 


are often developed from simplified models. It is important to note that real engine passages do not have perfect 


rectangular cross sections, but include corner fillet, ribs with fillet radii and special orientation. Therefore, this 


work provides detailed fluid flow and heat transfer data for a model of radial cooling geometry which possesses 


very realistic features. 
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Introduction 


Heat exchange is a critical issue in the development of 
gas turbine blades. Gas temperature at the first stage of 
the turbine is extremely high and can cause metal melting, 
thus, intensive cooling of the NGV and first blades is 
absolutely necessary. Usually a combination of internal 
convection cooling and film cooling is applied. They 
have to be designed for high efficiency, long life cycles 
and safe operation in spite of being exposed to high 
thermal loadings of the engine. 

The internal cooling systems of high pressure blades 
in modern aero engines are very complex and essentially 
rely on serpentine systems. The data available for design 
purposes tends to be in the form of correlations obtained 
from a large number of idealized geometries, e.g. [1, 2, 3] 
over restricted ranges of operation. The work described 
in the present paper is a part of the work undertaken 
within the EU ERICKA project [4]. The project con- 
cerned generation of test data on the heat transfer and 


Nomenclature 
X length of the channel (mm) 
Y height of the channel (mm) 


Z span of the channel (mm) 


pressure drop performance of ribbed cooling passages 
which, together with CFD and optimization tools, is in- 
tended to support the development of more effective 
cooling systems for rotating turbine components. 

Correlations for the heat transfer coefficient and pres- 
sure drop used in the design of radial cooling passages 
are often developed from data gained from simplified, 
flat surfaced models over specific conditions [5]. It is 
important to note that real engine passages do not have 
perfect rectangular cross sections, but include fillets and 
space oriented ribs also with fillet radii. Therefore, this 
work provides detailed fluid flow and heat transfer data 
for a model of radial cooling geometry which possesses 
very realistic features. The ribs are of a real shape. 

The main goal of this paper is to study the flow struc- 
ture and heat transfer, and also their correlations on all 
the four walls of a radial cooling passage model of a gas 
turbine blade. The flow structure and pressure drop were 
measured by means of the classical measurements tech- 
niques. The stagnation pressure and velocity were meas- 


PO inlet stagnation pressure (Pa) 
PP stagnation pressure (Pa) 
HTC, h heat transfer coefficient 


Nu Nusselt number 
Re Reynolds number 
d diameter (mm) 

T Temperature (K) 

k material constant 

u velocity (m/s) 
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f frequency (Hz) 
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ured in the channel outlet plane. The Transient Liquid 
Crystal method was used to investigate the heat transfer. 

The experimental setup presented here is able to pro- 
vide complementary data of the heat transfer and flow 
pattern in a cooled inner channel of a turbine blade. The 
transient heat transfer measurement technique using liq- 
uid crystal coatings is well known and has been success- 
ful in high resolution measurements on surfaces with 
complex shapes and geometries [6, 7, 8, 9]. The method 
allows obtaining detailed surface distributions of local 
heat transfer coefficients in the cooled channels, whereby 
it is possible to predict the result of the heat load under 
real-engine conditions. One of the advantages of this 
method is that there is no need to reach stable conditions 
in wall temperature, what is very important in the studies 
where duct geometry has a large surface and a highly 
curved and complex flow structure. 


Experimental Setup 


Heat transfer measurements 


To carry out research of heat transfer using the TLC 
technique [10, 11] it is necessary to have liquid crystal in 
the form of a sheet or sprayable microencapsulated water 
slurries placed on the inner wall of a test section. In this 
experiment a narrow band liquid crystal sheet (R35C1W) 
was used. The first step is calibration which is done by 
determination of the temperature at which the liquid 
crystal turns on a specific color. In this experiment, it was 
a green color. This is the way to determine the tempera- 
ture at which the green color intensity is the highest. The 
test section should be also equipped with thermocouples 
which measure the temperature along the flow and on the 
channel walls. In order to trigger the heat flux, it is ne- 
cessary to have a temperature difference between the gas 
and the channel walls. In this experiment the hot air was 
introduced into the channel triggering the heat transfer 
from the flowing gas to the surfaces of the channel walls. 
As the time was passing, the liquid crystals located on 
the wall surface were heated up and finally changed their 
color. Knowing the material properties of the channel 


LE leading edge 


TE trailing edge 
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walls, the local temperatures and also the time which 
passed from the start of the measurement (i.e. the mo- 
ment of introduction of hot air into the channel) to reach 
the wall temperature activating the green color, one can 
determine the heat transfer coefficient. The 3CCD cam- 
era recorded images during the measurement to deter- 
mine time ¢ at which the green intensity peak was 
achieved. The recorded frames were analyzed by means 
of our own automated MATLAB code. The software [12] 
for searching the time of the green intensity peak in every 
pixel was developed. Having the above information it 
was possible to determine the heat transfer coefficient A 
distribution on the wall. 

The liquid crystal calibration stand is shown in Fig. 1. 
The sheet with liquid crystals is placed on a copper plate, 
one of the edges of which is heated while the second 
edge is cooled using e.g. water at a constant temperature. 
The plate is placed in a housing which is intended to re- 
duce heat losses to the ambient. This kind of configura- 
tion provides an almost one-dimensional temperature 
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Fig. 1 Calibration of liquid crystal sheet. 
distribution along the plate. Thermocouples embedded in 
the plate measure the temperature at its surface. Then, 
combining the temperature distribution with green color 
distribution on the plate the temperature of green inten- 
sity peak can be determined. 

Fig. 2 shows the green color intensity versus time 
during the experiment for a chosen pixel. The tempera- 
ture at which the green color reaches the highest intensity 


is known from the calibration. Therefore, the surface 
temperature on the wall is determined exactly in the ex- 
periment. 
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Fig. 2 Green colour intensity history. 


Test section 


The view of the test section is shown in Fig. 3. This 
test section corresponds to an internal radial passage 
channel of a turbine blade. The lower and upper walls are 
ribbed, what corresponds to the pressure and suction side 
of the turbine blade inner channel. In the downstream 
view the left side wall corresponds to the leading edge 
direction, the right side wall corresponds to the trailing 
edge direction. Small ribs located on the lower and upper 
walls have a parallel setting to each other. The ribs were 
skewed by 45° in relation to the flow direction. The as- 
pect ratio of the channel is 1:4. The Reynolds number in 
the experiment based on hydraulic diameter was Re= 
63 000. 

The test section scheme together with the system of 
hot air preparation is presented in Fig. 4. The ambient air 
enters through the heater placed at the inlet of the settling 
chamber. Hot air is permanently evacuated from the set- 
tling chamber through the test section or the bypass sys- 
tem. While the bypass valve is open the air from the set- 
tling chamber flows to the bypass circuit. This arrange- 
ment is for establishing the mass flow and temperature, 
the same as during the measurement. When the bypass 
valve is closed, the main valve is open and hot air flows 
through the test section and then to the vacuum tank. The 
mass flow rate in the test section is controlled during the 
measurements by the orifice plate. For the investigation 
of heat transfer the liquid, crystal sheets were placed on 
the surface of the inner walls of the test section. 

The imaging system allows observing three walls, one 
sidewall and two ribbed walls simultaneously (see Fig.5). 
Two flat mirrors are mounted at 30 degrees and show 
upper and lower ribbed walls. Illumination is provided by 
four strip LED lamps. The white color of light and its 
precise orientation gives sufficient response of the green 
signal from the narrow band liquid crystal sheet. The 


imaging system comprises 3 CCD cameras which record 
images at a rate of 30 Hz and with a resolution of 1024 x 
768 pixels. 
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Fig.3 View of test channel. 
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Fig. 4 Measurement system scheme. 
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The local heat transfer coefficient is calculated in 
every pixel separately by solving the analytically unstea- 
dy, one dimensional conduction equation with the as- 
sumption of a semi-infinite wall for the whole of the test 
period, constant initial temperatures and a convective 
boundary condition. For this purpose a MATLAB code 
was specially created to calculate the heat transfer coeffi- 
cient h from the equation (1) developed by Clifford and 
Jones [10]. The MATLAB code is used to calculate the 
heat transfer coefficient at each pixel location by com- 
paring the experimental green intensity history with the 
gas temperature flowing through the test section. At the 
beginning the software calculates time ¢ after which there 
is the maximum intensity of the green color for every 
pixel separately. Then, using gas temperature Tgaw meas- 
ured by four thermocouples placed along the channel 
(linearly averaged) the local heat transfer coefficient is 
calculated. Ts). is the wall surface temperature at which 
liquid crystal adopts the green color. 7; is the initial wall 
temperature (constant during measurements) while k is a 
product of the channel material density, specific heat and 
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thermal conductivity. 
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Results and Discussion 


Pressure and velocity measurements 


The normalized stagnation pressure and velocity dis- 
tribution at the test section outlet plane (X=550, Fig. 3) is 
shown in Fig. 6. The results are shown according to the 
downstream view. The stagnation pressure was measured 
by means of a Pitot probe. The velocity was measured by 
means of a hot wire probe, hence, in this case the result 
shows the velocity magnitude. The X component (longi- 
tudinal velocity) is dominating in the channel flow. 

The pressure and velocity measurements show that a 
pair of main counterrotating vortices is created in the 
flow by the ribs. The dashed ellipses with an arrow show 
the sense of rotation of the induced vortices. The parallel 
setting of ribs in the channel creates two main vortices 
with the opposite sense of rotation, but the flow direction 
in these two vortices is the same in the channel centre, it 
is passing from one side wall to the other. In Fig. 6 one 
can also notice that regions of higher velocity correspond 
to regions of higher stagnation pressure. The highest flow 
velocity is in areas near the upper and lower walls, where 
the flow is forced by the skewed ribs, goes down over the 
rib toward the LE side wall. 


Heat transfer measurements 


Fig. 7 presents the heat transfer coefficient distribution 
on all four channel walls. The flow structure coming 
from aerodynamic measurements is placed on the right 


side. Arrows on the pressure distribution plane indicate 
the location in the space of particular walls. Fig. 7 shows 
that the HTC distribution is highly correlated with flow 
structure. 

Heat transfer strongly increases in areas where stagna- 
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tion pressure and velocity are high. The stagnation pres- 
sure is higher in the middle of the TE side wall than on 
the LE wall. This is reflected in the HTC distribution, in 
the middle of the duct, the HTC on the trailing edge wall 
(TE) is higher in relation to the leading edge wall. Corre- 
lation between pressure (velocity) and the heat transfer 
coefficient is also seen on the ribbed walls and in those 
channel corners where the flow has accelerated. The 
higher velocity in the left corners (at the leading edge 
wall) is observed, thus, in these regions one can see very 
high values of the heat transfer coefficient, in contrast to 
the opposite corner where HTC values are much lower. 
Generally one can say that higher heat transfer is in 
the region where shear stresses increase (Reynolds anal- 
ogy). This can be also illustrated zooming the region be- 
tween the ribs (Fig. 8). The ribs themselves are not cov- 
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Fig. 6 Measurement plane downstream of ribbed section, a) 
normalized stagnation pressure, b) velocity magnitude. 
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Fig. 7 Heat transfer coefficient distribution on all channel walls. 


ered by LCs. The shear stresses increase and the heat 
transfer coefficient also significantly increases just 


downstream of the rib where the flow driven by the vor- 
tex goes down. The heat transfer radically decreases in 
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the region before the ribs where the flow separates. 

Watching the HTC distribution on the leading edge 
side wall (Fig. 7) one can also observe the phenomenon 
of vortices development along the channel. At the inlet 
where the ribs start to influence the flow the height of 
two main vortices is small, then, along the channel one 
can see that they gradually enlarge. Finally at the outlet 
the vortices occupy the whole cross-section of the channel. 

The averaged normalized Nusselt number (Nu/Nu0) 
distributions along the channel for side and ribbed walls 
are shown in Fig. 9. For side walls the averaging was 
made calculating the mean value, vertically line by line 
along the channel. For the ribbed walls Nusselt number 
was calculated averaging the surface between the neigh- 
boring ribs. 


Fig. 8 Heat transfer coefficient distribution on ribbed wall. 


Fig. 9 shows that the heat exchange increases on the 
side walls along the channel from the inlet to about 80% 
of the duct length. At the outlet the average Nusselt 
number values start to decrease. The parallel setting of 
ribs causes the difference in Nu distribution between LE 
and TE walls. On the LE wall, where the direction of 
vortices is out from the wall (see Fig. 7), lower average 
Nusselt number values in relation to TE wall are ob- 
served. 
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Fig. 9 Averaged normalized Nusselt number distribution on 
channel walls. 


Figure 9 also shows the heat exchange along the 
ribbed walls on the pressure side and suction side. In this 
figure one can see that heat transfer increases quickly 
from the very beginning, what is due to vortices devel- 
opment in the channel. When the vortices have developed 
the average Nu stabilizes, further downstream of the 
channel there is a small tendency in the average heat ex- 
change to be decreasing in the direction of the duct outlet. 
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The vortices generated by the ribs increase significantly 
the heat exchange on the walls. The averaged Nusselt 
number values on the ribbed walls are approximately 50% 
higher in relation to the side walls. At the outlet of the 
channel one can also observe the difference in the heat 
transfer between the suction (SS) and pressure (PS) sides. 
The reason of this phenomenon is a non-uniform velocity 
distribution in the channel. In Fig. 6b one can see that 
velocity at the upper wall (SS) is a little higher than at the 
lower wall (PS). Lower velocity means lower shear 
stresses which in consequence lead to a decrease in the 
heat transfer coefficient. On this example one can also 
conclude that heat transfer can change due to the local 
small asymmetry of the flow field. 


Velocity fluctuation measurements 


The result of the velocity fluctuations measurement by 
HTA (Hot Wire Anemometry) is shown in Fig. 10. The 
measurement plane was located downstream of the 
ribbed section at X=550, in the same plane adequate for 
the velocity and stagnation pressure measurements. 
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Fig.10 Normalized velocity fluctuation. 


The results in Fig. 10 show that the lowest level of the 
velocity fluctuations in the flow can be observed in the 
acceleration regions near the bottom and upper wall of 
the channel. The highest velocity fluctuations are observed 
near the side walls in areas where the flow decelerates. 
Generally, we can say that a high level of turbulence can 
be observed in the channel with the aspect ratio 1:4, 
which is relatively narrow. The average turbulence inten- 
sity is at a level of about 10%. A high level of turbulence 
on the one side supports the heat exchange, but slightly 
increases the losses in the channel flow on the other side. 


Conclusions 


The following conclusions can be drawn from the 

above described investigations: 

e A model of an internal cooling channel of a gas tur- 
bine blade was manufactured and equipped with 
imaging and necessary systems for heat transfer and 
aerodynamic measurements; 

e The method of heat exchange measurements in a 
cooling channel of a gas turbine blade was based on 
the transient liquid crystal technique. The transient 
method is not expensive in implementation and its 
biggest advantage is that there is no need to achieve 
stable conditions in wall temperature; 

e A comparison of heat exchange with aerodynamic 
measurement (flow pattern) in the ribbed channel 
was made; 

e Ribs on the suction and pressure side walls create a 
specific flow structure with two dominant vortices 
along the channel with the opposite sense of rota- 
tion; 

e A high level of turbulence (average turbulence in- 
tensity of about 10%) in the ribbed and narrow 
channels was observed; 

e The heat transfer coefficient distribution on the 
walls is highly correlated with the vortex pattern 
and the stagnation pressure distribution in the 
channel; 

e The heat transfer strongly increases in areas where 
local velocity (stagnation pressure) is higher. On the 
side walls these are areas where the flow is directed 
towards the wall. On the ribbed walls these are 
areas just downstream of the rib and this side of the 
corner where the flow in the vortex is upcoming to 
that corner. The heat transfer decreases considerably 
upstream of the next rib and the opposite corner; 

e In the radial cooled passage channel of the turbine 
blade the differences in the heat transfer between 
the suction side and pressure side can appear due to 
the local small asymmetry of the flow field. 
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